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Abstract. A gaseous Bose-Einstein condensate (BEC) offers an ideal testing ground for studying 
symmetry breaking, because a trapped BEC system is in a mesoscopic regime, and situations exist 
under which symmetry breaking may or may not occur. Investigating this problem can explain why 
mean-field theories have been so successful in elucidating gaseous BEC systems and when many- 
body effects play a significant role. We substantiate these ideas in four distinct situations: namely, 
soliton formation in attractive BECs, vortex nucleation in rotating BECs, spontaneous magnetization 
in spinor BECs, and spin texture formation in dipolar BECs. 
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l> : INTRODUCTION 



0\ 

Whereas the ground state of a microscopic system possesses the exact symmetry of the 
Hamiltonian, the ground state of a macroscopic system seldom, if ever, represents the 
full symmetry of the Hamiltonian. This phenomenon, known as spontaneous symmetry 
breaking, is the key concept for understanding macroscopic phenomena. The purpose of 
this study is to demonstrate that a gaseous Bose-Einstein condensate (BEC) offers a new 
paradigm for studying this problem. 

Many insights into the nature of superfluidity have been gained through the study 
of superfluid helium systems. It is therefore instructive to compare the properties of 
atomic-gas BECs with those of superfluid helium. 

One major difference is the kinetics. The collision time for liquid helium is much 
^ shorter than the inverse collective mode frequency. This implies that the local thermo- 

dynamic equilibrium is ensured and that the physics can be understood by conservation 
laws and hydrodynamics. In the case of an atomic-gas BEC, the collision time is of the 
order of the inverse collective mode frequency. This implies that the local thermody- 
namic equilibrium cannot always be achieved and that once the system is driven out of 
equilibrium, the ensuing non-equilibrium relaxation and kinetics are essential for under- 
standing such phenomena as BEC phase transition and vortex nucleation. 

Another major difference is symmetry breaking. In the case of bulk liquid helium, 
the thermodynamic limit is always achieved, and spontaneous symmetry breaking of the 
relative gauge occurs, resulting in the emergence of the mean field. In the case of an 
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FIGURE 1. Stationary solutions of Eq. Q for several values of y. When y < 1, the ground state is a 
uniform condensate, and when / exceeds 1, the ground state develops into a soliton. 



atomic-gas BEC, the system is in the mesoscopic regime and the thermodynamic limit 
is, in general, not achieved. Therefore, whether or not the symmetry breaking occurs 
depends on the situation of the system, and if it does, its dynamics should be observable 
due to the aforementioned long collision time. 

In the present paper we substantiate these ideas in four distinct cases: namely, soliton 
formation in quasi one-dimensional attractive BECs, where the translational symmetry 
is broken; vortex nucleation in rotating BECs, where the axisymmetry is broken; spon- 
taneous magnetization in spinor BECs, where the rotational and chiral symmetries are 
broken; and the spontaneous formation of spin textures in spinor dipolar BECs, where 
we can expect the Einstein-de Haas effect and the ground-state mass flow. 



SOLITON FORMATION IN A QUASI-1D ATTRACTIVE BEC 

We begin by discussing symmetry breaking in a quasi-lD attractive BEC HI]. Suppose 
that an attractive BEC is confined in a quasi- ID torus. On a mean-field level, the 
properties of the system can be described by the Gross-Pitaevskii (GP) equation 

(-^-^ol 2 )^^, (1) 

where x Vq{0) is the ground-state wave function and 7 is the dimensionless strength of 
interaction. Figure 1 shows ground-state wave function as a function of the angular 
coordinate 6 for various strengths of interaction 7. 

When dimensionless strength of interaction 7 is smaller than 1, the ground-state den- 
sity is uniform. However, when it exceeds 1, the translational symmetry is spontaneously 
broken and a bright soliton is formed as shown in Fig. 1 . Thus the mean-field theory pre- 
dicts a second-order quantum phase transition at 7 = 1 . 

Figure 2(a) shows the Bogoliubov spectrum before and after breaking the translational 
symmetry |y]. The zero-energy mode emerges in the soliton regime, this being the 
Goldstone mode associated with the breaking of translational symmetry. Figure 2(b) 
shows the corresponding many-body spectrum [2]. The many -body spectrum in the 
uniform BEC regime is similar to that of the Bogoliubov spectrum. However, a dramatic 
change in the landscape of the energy spectrum occurs in the bright soliton regime. 




In particular, a quasi-degenerate spectrum appears above the ground state in the bright 
soliton regime. This quasi-degeneracy is a signature of the breaking of the translational 
symmetry that generates a bright soliton. It is remarkable that many-body physics can 
automatically generate such a symmetry-breaking-inducing quasi-degenerate spectrum 
which is absent at the mean- field level. 



AXISYMMETRY BREAKING IN A ROTATING BEC 

Let us next discuss symmetry breaking in vortex nucleation. When the rotational fre- 
quency of the container exceeds a certain critical value, the BEC begins to rotate and a 
vortex enters the system. The axisymmetry of the system is spontaneously broken in this 
vortex nucleation process. 

Figure 3 shows the many-body excitation spectrum of a rotating BEC as a function of 
the angular momentum L of the system with N particles, in which the energy is measured 
from the lowest energy state for each angular momentum J3]. The ground state is quasi- 
degenerate when the rotation frequency is low. This quasi-degeneracy may be regarded 
as a precursor for spontaneous symmetry breaking due to vortex nucleation. We note 
that the energy difference between these quasi-degenerate levels is of the order of l/N. 
This quasi-degeneracy is therefore solely of many-body nature and should vanish at 
the thermodynamic limit. In fact, we can show that the Goldstone mode associated 
with the axisymmetry breaking is the lowest-lying envelope of this quasi-degenerate 
spectrum yD. 

It is striking that a many-body quasi-degenerate state is spontaneously generated when 
the continuous symmetry of the system is about to break and then removed when its role 
is over, i.e., after the symmetry has been broken. This mechanism for symmetry breaking 
appears universal, regardless of the details of the underlying physics, and manifests itself 
only in the mesoscopic regime. 
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FIGURE 3. Many-body energy spectrum of a rotating BEC as a function of angular momentum L of 
a system with N particles. The energy is measured from the lowest-energy state for each L, and excited 
states that involve center-of-mass motion are not shown. The envelope of the lowest-lying excited states 
that is labeled G represents the Nambu-Goldstone mode associated with axisymmetry breaking. As L 
increases, the envelope develops an energy gap of the order of N, indicating mass acquisition by the 
Nambu-Goldstone mode [3]. 



SPINOR BEC: SPIN TEXTURE AND CHIRAL SYMMETRY 

BREAKING 

A long-standing question with magnetism is how the spontaneous magnetization of a 
ferromagnet can occur in an isolated system in which the total spin angular momentum is 
conserved. One possible solution is that all spins align in the same direction and that the 
system is in a quantum-mechanical superposition state over all directions. We propose 
here a different scenario, in that the system develops local magnetic domains of various 
types that depend on the nature of the interaction, conservation laws, and the geometry 
of the trapping potential. 

As an example, let us consider a spin-1 ferromagnetic BEC in a cigar-shaped trap. We 
assume that almost all the atoms are initially prepared in the magnetic sublevel m = 0. 
As time progresses, the m = population is converted into the m = ±1 modes due to the 
spin-exchange interaction, and the system develops spin textures. 

Figure 4 shows the time evolution of the mean spin vector on the trap axis 0]. The 
length of the spin vector is zero everywhere because of the initial condition. As time 
progresses, staggered magnetic domains develop due to dynamical instability. After 
more time has elapsed, helical structures are spontaneously formed [Fig. 4(d)] because 
the kinetic energy is decreased by the formation of a helix. 

The ferromagnetic BEC yields yet another surprise: chiral symmetry breaking. Sup- 
pose that all atoms are prepared again in the m = state in a pancake- shaped trap. Fig- 
ure 5 shows the Bogoliubov spectrum as a function of dimensionless strength (oc N) 
of the spin-exchange interaction @]. In the region of g^ D < —3.9, the modes with orbital 
angular momentum £ = ±1 have imaginary parts, and they are therefore dynamically un- 



(a) (b) (c) (d) 




co ± t=440 480 520 630 



FIGURE 4. Mean spin vectors along the z axis seen from the — y direction. The trap frequencies are 
Ox — 2?r x 245 Hz and co z — 2n x 4 Hz. A mag netic field of 54 mG is applied in the z direction. The 
initial state is given by = 0, y/b = VT — 10~ 4 y/ g , and = 0.01 y/ g , where y/ g is the m = — 1 ground 
state. The size of each strip is 600 x 16 in units of (h/mco^) 1 / 2 ~ 0.69 /im. The spin vector is displayed 
from z = jum to z = 54 jim in (a), (b) and (c), and from z = 35 fjm to z = 89 /im in (d). All the vectors 
are shown in the frame rotating about the z axis at the Larmor frequency. 



stable and grow exponentially. When a Rb BEC is prepared in this region, the m = 
atoms are transfered into the m = ±1 states due to the dynamical instability, and they 
obtain the orbital angular momentum. The angular momentum conservation implies that 
the m = 1 and m = — 1 components must have opposite sign of orbital angular momen- 
tum. There are two possibilities: the m = 1 component can have either orbital angular 
momentum of£=lor£ = — 1; correspondingly, the m = — 1 component can have either 
— 1 or 1 orbital angular momentum. These two possibilities are degenerate, this degen- 
eracy being a statement of the chiral symmetry. However, since a chirally symmetric 
state has higher energy than a chiral- symmetry -broken states [5|, the chiral symmetry is 
dynamically broken and each spin component begins to rotate spontaneously. 

Figure 5(b) shows the time development of the orbital angular momentum of the 
m = — 1 component o]. This remains zero for a certain latency period and then acquires 
a non-zero value due to the chiral symmetry breaking. Figures 5(c) and (d) show the spin 
textures for chirally symmetric and broken symmetry states o]. The chirally symmetric 
state has a domain wall at the cost of the ferromagnetic energy, while the chiral- 
symmetry-broken state circumvents this energy cost by developing topological spin 
textures. This represents the underlying physics of this chiral symmetry breaking. 

Our predictions have recently been observed by the Berkeley group |6]. They carried 



FIGURE 5. (a) Real and imaginary parts of lowest Bogoliubov energies gW for I = and ±1, where 
label I indicates that the m = ±1 components of the eigenfunction are proportional to e ±l ^ with (r, <j>) 
being the cylindrical coordinates in two-dimensional system. The two energy levels of e^ 1 ) are degenerate 
due to the axisymmetry of the system. We have taken the parameters of spin-1 87 Rb atoms, where spin- 
independent interaction strength gQ D is related to spin-dependent strength g^ D by gQ D = — 216.1g] D . (b) 
Time development of the fraction w_j and orbital angular momentum per particle L_i in the m = — 1 
component with (7 = 0.03, dashed) and without (7 = 0, solid) dissipation. The interaction strengths are 
gf = 2200 and g\ D = -10.18. The initial state is given by y/ = XfTg, = l0' A r{e^ +1.0001e^)v^, 
and y/\ = 0, where is the m = Q mean-field ground state, (c) and (d) show snapshots of the spin textures. 
The size of the frame is 12 x 12 in units of a^ . 



out experiments subject to the same initial conditions, i.e., all the atoms were initially 
prepared in the m = state. As we predicted, the system remained unmagnetized during 
a certain latency period, before spontaneously developing magnetization. They also 
observed a polar-core vortex corresponding to our chiral- symmetry -broken state. 

DIPOLAR BEC: EINSTEIN-DE HAAS EFFECT AND 
GROUND-STATE MASS FLOW 

Let us finally discuss a dipolar BEC by focusing on the spin degree of freedom. The 
magnetic dipole-dipole interaction is a tensor force which causes spin-orbit coupling, 
and only the total, spin plus orbital, angular momentum is conserved. Therefore, the 
dipolar interaction transfers angular momentum between the spin and orbit, i.e., the 
Einstein-de Haas effect occurs. 

The Stuttgart group created a spin-polarized 52 Cr BEC under a relatively strong mag- 
netic field fZD. What would happen to this BEC if we reduced the magnetic field to zero? 
We could expect spin relaxation to occur, and because of the conservation of total angu- 
lar momentum, BEC would start to rotate spontaneously. This is the Einstein-de Haas 
effect To check this idea, we performed numerical simulations by solving the 

seven-component non-local GP equation in three dimensions. We find that the BEC be- 
gins to rotate spontaneously to compensate for a decrease in the spin angular momentum, 
as shown in Fig. 6(a). We also find that atomic spins undergo Larmor precession around 
the local dipole field, and develop topological spin textures as shown in Fig. 6(b) and 
(c). In this case, the dipole field is directed inward in the upper hemisphere and outward 
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FIGURE 6. (a) Time development of the spin and orbital angular momenta in a spherical trap in the 
absence of a magnetic field. The initial state is spin polarized at m = —3. (b) and (c) show snapshots of 
the spin textures on (b) the z — 2 fim plane and (c) the z — —2 fim plane at t — 0.4 ms, where the Thomas- 
Fermi radius is 3.5 jtim. The length of the arrow represents the magnitude of the spin vector projected on 
the x—y plane. 



in the lower hemisphere, and we therefore have spin textures with opposite directions. 

The dipole-dipole interaction is also expected to yield ground- state spin texture in a 
ferromagnetic BEC, as in the case of the domain structure in a solid-state ferromagnet. 
The Einstein-de Haas effect is also known to occur in a ferromagnet. Then, what is truly 
new in spinor dipolar BECs? 

The unique feature of the spinor dipolar BEC which is absent from a solid-state 
ferromagnet is the spin-gauge symmetry which relates the spin texture to a mass current. 
The fundamental query is whether or not a spinor dipolar BEC can exhibit a spontaneous 
mass current in the ground state. The answer is summarized in the phase diagram shown 
in Fig. 7(a) lllflll . in which i?TF is the Thomas-Fermi radius, <^dd is the dipole healing 
length, and <^ sp is the spin healing length. 

We have found three phases. When the system size is small, the flower phase is sta- 
bilized in which spins are almost polarized and slightly flare out depending on RTp/^dd- 
When the dipolar interaction becomes stronger, or equivalently the system size becomes 
larger, the spin vectors tilt into the either (j) or — direction, where (j) is the azimuthal 
angle. This spin texture has chirality and we call this phase the chiral spin- vortex phase. 
In this phase, BEC has a substantial net mass current which monotonically increases as a 
function of Rtf/^^ leading to increase of the kinetic energy. Therefore, when RTp/^dd 
becomes larger, the system goes to the next phase: the polar-core vortex phase, the spin 
texture in this phase forming a vortex with no net circulation. 

All three phases can be realized in the spin-1 87 Rb BEC by varying the trap frequency. 
Figure 7(b) shows the orbital angular momentum of the ground-state spin-1 87 Rb BEC 
as a function of the trap frequency and the number of atoms. We find particularly in the 
chiral spin-vortex state, that the orbital angular momentum is increased by up to 40 % 
of the full value of the singly quantized circulation. 
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FIGURE 7. (a) Phase diagram of a spin-1 ferromagnetic dipolar BEC. The solid curve shows the 
boundary between the states with total angular momentum 7 = and / = 1, and the broken line divides 
the phases with and without spin chirality. The schematic diagram shown in each region represents the 
spin configuration (black arrows) and mass flow (gray arrows). The spin-1 87 Rb BEC traces the dotted 
line where ^ S p/<^dd = 0.30. (b) Orbital angular momentum as a function of trap frequency co and number 
of atoms N for a spin-1 87 Rb BEC along the dotted line in Fig. 7(a). In the case of a BEC with a million 
atoms, the chiral spin-vortex phase exists in the region of 2n x 70 Hz < 0) < 2n x 630 Hz. 



SUMMARY 

A gaseous BEC is a mesoscopic system. Studying this system enables us to learn in detail 
when and how symmetry breaking occurs, and to identify the dynamics of symmetry 
breaking because of the long collision time. Various types of symmetry breaking are 
experimentally accessible. In the case of soliton formation, the translational symmetry 
is spontaneously broken; in the case of vortex nucleation, the axisymmetry is broken; in 
the case of a spinor BEC, the rotational and chiral symmetries are broken and various 
spin textures are spontaneously generated due to competition between the conservation 
law, interaction, and geometry of the system; and in the case of a dipolar BEC, we 
expect the Einstein-de Haas effect and substantial ground-state circulation in the chiral 
spin- vortex phase. 
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